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Appendix  B 


I ntroduct i on 

This  Appendix  contains  the  graphical  results  of  the 
Monte  Carlo  analysis  of  this  study.  The  plots  will  be  pre- 
sented in  sets.  All  sets  will  include  the  dynamic  state 
error  plots  for: 

V*  - x-velocity  for  the  missile  in  the  inertial  frame 
mx 

0T  - Line -of-sight  angle  as  seen  by  the  ownship 
R - Range 

R - Range  rate  (relative  closing  velocity  where,  by 
convention,  R = -V 

c 

- Missile’s  developed  lateral  acceleration. 

And  when  estimated,  the  parameter  error  plots  for : 
n - Proportional  navigation  constant 
Tjp.  - Time  constant  of  the  first  order  lag 
M/S  - Ratio  of  mass  to  cross  sectional  area 
will  be  included  in  the  set.  In  addition  to  the  error  plots, 
a tuning  plot  (covariance  matching)  will  be  included  for  each 
state  and  parameter.  The  smooth  curve  on  the  tuning  plot 
represents  the  square  root  of  the  filter-calculated  covar- 
iance. The  "apparently  noisy"  curve  represents  the  square 
root  of  the  variance  of  the  error  calculated  over  the  twenty 
simulations.  Also  included,  in  various  sets,  are  covariance 
convergence  plots  (described  in  Chapter  IV).  The  box  in  the 
upper  right  hand  cornel1  of  each  plot  indicates  whether  the 
plot  is  a mean  error  plot  or  a covariance  plot.  The  order 
of  the  missile  a:,  well  as  the  data  set  number  are  also  in- 

1 

# J 

A 


eluded.  The  data  set  number  indicates  the  tuning  parameters 


used  for  a particular  run.  The  filter  initial  estimates,  the 


tuning  parameters  and  a brief  description  of  the  purpose  of 


each  set  v/ill  be  included  at  the  beginning  of  each  set.  It 


should  be  noted  that  the  high-g  scenario  was  used  for  all 


cases  except  the  one  annotated  low-g. 


Organization 


The  sets  v/ill  be  arranged  in  the  following  order: 


Figure 


Subject 


Zero  order  missile  filter 


.10-21 


Third  order  missile  filter 


22-31 


Benchmark 


32-41 


Fourth  order  missile  filter  (A/P  at  0 secj 


42-51 


Fourth  order  missile  filter  (A/P  at  3 sec) 


52-61 


62-71 


Fourth  oi'der  missile  filter  (A/P  at  5 sec) 
Fourth  order  missile  filter  (A/P  at  0 sec  - 


72-81 


complete  linearization  of  f) 

Fourth  order  missile  filter  (A/P  at  5 sec  - 


using  fourth  order  Runge-Kutta  integration 


82-91 


package ) 

Sensitivity  analysis  (n  = 6.) 


92-101 


Sensitivity  analysis  (n  = 3. ) 


J 02-111 


Sensitivity  analysis  (r2  = *8) 


112-121 


Sensitivity  analysis  (r2  = -1) 


122-131 


Sensitivity  analysis  (M  = 2.) 


132-141 


Sensitivity  analysis  (M  = 8.) 


142-151 


First  order  missile  filter  (rf  set  equal  to  .89) 


2 


152-164 

n estimation  - n initialized  at  3- 

165-178 

n estimation  - n initialized  at  6. 

179-190 

Tf  estimation  - initialized  at  1.5 

191-202 

estimation  - initialized  at  .3 

203-216 

M/S  estimation  - M/S  initialized  at  45. 

217-230 

M/S  estimation  - M/S  initialized  at  15- 

231-244 

n and  estimation  (high-g  scenario) 

245-258 

n and  estimation  (low-g  scenario) 

259-276 

n,  t^,  and  M/s  estimation 

277-294 

n,  and  M/S  estimation  with  the  dyna: 

states  initialized  with  some  error. 
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Zero  Order  Missile  Filter 


The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  arc 


VT  (0) 
mx'  ' 

= 1225-7  fps 

0(0) 

= 4.  3633^5  radians 

R(0) 

= 10000.  feet 

R(0) 

= -2122.  fps 

a.  ( 0 ) 

j J 

= 0. 

n(  0) 

-4.5 

Tf(0) 

= n/a 

M/S(0)  = 29.197  siugs/ft2 


These  plots  wore  generated  to  demonstrate  the  per- 
formance of  this  filter  with  the  software  changes  noted  5 r 
Appendix  A.  The  parameters  n and  M/S  were  not  estimated, 
but  were  set  to  their  correct  values. 
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Fig.  4.  ANCLE  MEASUREMENT  SICMAS  ZERO  ORDER 
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Fig.  7.  RflNCE  RRTE  ZERO  ORDER  MISSILE 
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Third  Order  Missile  Filter 


The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 


yL<°> 

= 1225.7  fps 

9(0) 

= 4.363345  radians 

R(0) 

= 10000.  feet 

R(0) 

= -2122.  fps 

aL(0) 

= 0. 

n(0) 

= 4.5 

rf(0) 

= N/A 

M/S(0) 

= 29.197  Slugs/ft2 

R = 


T.E-5  0.  0. 
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This  design  was  used  as  an  attempt  in  establishing  a 
benchmark.  The  filter  model  included  the  exact  structure  of 
the  autopilot  (found  by  assuming  constant  coefficients  as  de- 
scribed in  Chapter  III).  Chapter  III  lists  six  sets  of  co- 
efficients found  at  various  times  over  the  high-g  scenario. 
The  plots  that  follow  were  generated  using  the  coefficients 
for  t=0.  This  set  of  coefficients  reduced  the  initial  tran- 
sients in  the  filter.  This  can  be  seen  by  comparing  the 
plots  of  the  fourth  order  missile  filter  which  use  three 
sets  of  autopilot  transfer  function  coefficients  to  demon- 
strate this  point. 
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Fig.  14.  RNOLE  MEASUREMENT  SI  OMRS  THIRD  ORDER 


Fig.  16.  RRNGE  SlGMHo  THIRD  ORDER 
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Pig.  19.  LRT  RCCLLERRT I ON  THIRD  ORDER  MISSILE 
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Fig.  ?0.  LRT  RCCELERRTIQN  SI  OMRS  THIRD  ORDER 
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r : 

Benchmark 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 

V1  (0)  = 1225.7  fps 

mx 

9(0)  = 4 . 363345  radians 

R( 0 ) = 10000.  feet 

R(0)  = -2122.  fps 

a^( 0)  - 0. 

n(0)  = 4.5 

Tf(0)  = n/A 

M/S(0)  = 29.197  slugs/ft2 
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A 

complete  des 

cription  of 

the  £d 

hoc  de 

sign, 

used  in 

generating  the  following 

plots, 

can  be 

found 

in  Ch 

apter  III 

Basically,  the  true 

commanded  accelera 

tion  of  the 

missile 

was  assumed  known  perfectly  and 

passed 

direc 

tly  to 

the  auto 

pi]ot  in  the  filter 

model 

. The 

autopilot  transfer 

function 

used  the  coefficients  for  t=0  (refer  to  Chapter  III).  These 
plots  were  considered  as  benchmarks  for  the  five  dynamic 
states  of  the  fundamental  filter  design. 
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X-VELOCITY 


Fig.  22.  X-VELOCITY  THIRD  ORDER  MISSILE 
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Fig.  26.  RANGE  THIRD  ORDER  MISSILE 
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Fourth  Order  Missile  Filter  (A/p  at  0 sec ) 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 

V1  (0)  = 1225.7  fps 
rax 

0(0)  = ^.3633^5  radians 

T 

R( 0)  = 10000.  feet 

R( 0)  = -2122.  fps 

a^o)  = 0. 

n( 0)  =4.5 

Tf(0)  = N/A 

M/S ( 0 ) = 29.197  Slugs/ft2 


These  plots  were  generated  by  the  fourth  order  filter 
with  the  coefficients  of  the  autopilot  transfer  function  de- 
termined at  t=0.  These  plots  were  compared  to  those  results 
of  the  fourth  order  filter  using  autopilot  coefficients  for 
t~3  and  t=5.  The  intent  was  to  determine  if  there  was  any 
distinguishable  difference  in  filter  performance  for  the 
various  sets  of  autopilot  coefficients.  A complete  descrip- 
tion of  the  filter  autopilot  and  its  coefficients  can  be 
found  in  Chapter  111. 
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Fig.  32.  X-VELOCITY  FOURTH  ORDER  MISSILE 
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Fig.  34.  angle  measurement  fourth  order  missile 


45 


5.35 


Fig.  35-  ANGLE  MEASUREMENT  SIGf.AS  FOURTH  ORDER 
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Fig.  39-  RANGE  RATE  SIGMAS  FOURTH  ORDER 
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Fig.  40.  LHT  RCCELERRT10N  FOURTH  ORDER  MISSILE 
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Fig.  4l.  LRT  RCCELERRT10N  SICNRS  FOURTH  ORDER 
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Fourth  Order  Missile  Filter  ( A/P  at  sec ) 


The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 

Vmx(0)  = 1225,7  fpS 
9^(0)  = 4.3633^5  radians 

R ( 0 ) = 10000.  feet 

R ( 0 ) = -2122.  fps 

ar  ( 0 ) = 0 . 

n(0)  = 4.5 

Tf(0)  = N/A 

M/S( 0)  =29.197  Slugs/ft2 
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These 

plots  v;e  re  generated 

by  the  fourth  order 

f il ter 

with  the  cc 

efficients  of 

the  autopilot  trans 

fer  function 

de- 

termined  at 

t=3.  These 

plots  were  compared 

to  thos 

e resul ts 

of  the 

fourth  order  filter  usin, 

g autopilot  coeffici 

ents  i 

"or 

t~0  and  t=5 

. The  intent 

was  to 

deternine  if 

there 

was  any 

distinguishable  difference  in  filter  performance  for  the 
various  sets  of  autopilot  coefficients.  A complete  descrip- 
tion of  the  filter  autopilot  and  its  coefficients  can  be 
found  in  Chapter  III. 
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Fig.  43.  X-VELDCITY  SIGNRS  FOURTH  ORDER 
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Fig.  48.  RfiNCE  RRTE  FOURTH  ORDER  MISSILE 
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Fig.  49.  RANGE  RATE  SIGMAS  FOURTH  ORDER 
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Fig.  50.  LRT  RCCELERRT1 ON  FOURTH  ORDER  MISSILE 
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Fourth  Order  Missile  Filter  ( A/P  at  sec ) 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 


V1  (0)  = 1225.7  fps 
mx 

0(0)  = 4.363345  radians 

T 

R(0)  = 10000.  feet 

R ( 0 ) = -2122.  fps 

aL(0)  = 0. 

n(0)  = 4.5 

Tf(0)  = N/A 

M/S( 0)  = 29.197  slugs/ft2 
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0; 

These  plots  were  generated  by  the  fourth  order  filter 
with  the  coefficients  of  the  autopilot  transfer  function  de- 
termined at  t=5.  These  plots  were  compared  to  those  results 
of  the  fourth  order  filter  using  autopilot  coefficients  for 
t=0  and  t=3-  The  intent  was  to  determine  if  there  was  any 
distinguishable  difference  in  filter  performance  for  the 
various  sets  of  autopilot  coefficients.  A complete  descrip- 
tion of  the  filter  autopilot  and  its  coefficients  can  be 
found  in  Chapter  III. 
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Fig.  53.  X-VELOCITY  SIUMRS  FOURTH  ORDER 
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Fig.  54.  RNDLE  MERSUREMENT  FOURTH  ORDER  MISSILE 
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Fig.  55.  ANGLE  MEASUREMENT  SIGMAS  FOURTH  ORDER 
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Fig.  56.  RANGE  FOURTH  ORDER  MISSILE 
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These  plots  were  generated  to  justify  the  assumptions, 
that  9rn,  V , and  Y could  he  considered  constant  over  the 
measurement  period  of  0.02  seconds.  The  fourth  order  filter 
was  used  with  the  coefficients  of  the  autopilot  chosen  for 
t=0.  The  F-matrix  within  the  filter  contains  a complete 
linearization  of  the  f-vector.  A comparison  of  these  p^ots 
with  those  for  "Fourth  Order  Missile  Filter,  A/P  at  t-0" 
(Figures  32  through  4l),  can  be  made  to  determine  the  dif- 
ference v/hen  the  above  mentioned  variables  are  considered 
constant  over  the  measurement  period. 
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Fig.  62.  X-VEL0C2TY  FOURTH  ORDER  MISSILE 
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Fig.  64.  RNOLE  MEASUREMENT  FOURTH  ORDER  MISSILE 
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Fig.  65.  ANGLE  MEASUREMENT  SIGMAS  FOURTH  ORDER 
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Fourth  Order  Fissile  Fii  ter  (using  Rur.ge-Kutta  integration) 
The  initial  state  estimates  and  the  tuning  parameters 

for  this  case  are 


ymx<°> 

= 1225.7  fps 

0(0) 

= 4.363345  radians 

R(0) 

= 10000.  feet 

R(0) 

= -2122.  fps 

aL(0) 

= 0. 

n(  0) 

= 4.5 

Tf(0) 

= N/A 

M/S(0)  = 29.197  slugs/ft2 
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These  plots  were  generated  to  determine  if  the  filter 
performance-  was  degraded  by  using  a first  order  Euler  inte- 
gration routine.  The  results  shown  here  used  a Runge-Kutta 
routine  and  can  be  compared  to  Figures  52  through  6l,  which 
were  generated  with  the  Euler  routine.  The  coefficients  of 
the  autopilot  were  chosen  for  t=5.  since  these  coefficients 
caused  the  greatest  transients  in  the  dynamic  states  of  the 
filter  (compared  to  t=0  and  t=3.  as  described  in  Chapter  III). 
It  was  decided  that  a favorable  comparison  for  this  case 
would  be  the  best  test  for  the  Euler  routine. 
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Fig.  77.  RHNGE  SIGMRS  FOURTH  ORDER 
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Fig.  80.  LRT  flCCElERflT ION  FOURTH  ORDER  MISSILE 
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Sensitivity  Analysis  (n  = ) 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 
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This  set  of  plots  was  generated  by  setting  the  propor- 
tional navigation  constant,  n,  to  6 in  the  truth  model.  The 
fourth  order  filter  was  used  with  n set  to  4.5  in  the  filter. 
Only  the  dynamic  states  of  the  missile  model  were  estimated 
by  the  filter. 
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Sensitivity  Analysis  (n  ~ 3.  ) 


The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 
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This  set  of  plots  was  generated  by  setting  the  propor- 
tional navigation  constant,  n,  to  3 in  the  truth  model.  The 
fourth  order  filter  was  used  v/ith  n set  to  4.5  in  the  filter. 
Only  the  dynamic  states  of  the  missile  model  were  estimated 
by  the  filter. 
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Fig.  93.  X-VEL0C1 TY  SIDMRS  FOURTH  OROER 
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Fig.  100.  LRT  RCCELERRT ION  FOURTH  ORDER  MISSILE 
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Fig.  101.  LOT  ACCELERATION  SIGMAS  FOURTH  ORDER 
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Sensitivity  Analysis  ( = . 3 ) 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 
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This  set  of  plots  was  generated  by  setting  the  tine 
constant  of  the  guidance  system,  r2  to  seconds  in  the 
truth  model.  r2  in  the  filter  was  set  to  0.3  seconds.  The 
fourth  order  filter  was  used  and  only  the  dynamic  states  of 
the  missile  model  were  estimated. 
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Fig.  108.  RRNGE  RRTE  FOURTH  ORDER  MISSILE 
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Fig.  110.  LRT  fiCCELERRTIDN  FOURTH  ORDER  MISSILE 
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Sensitivity  Analysis  (r  - .1) 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 
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This  set  of  plots  was  generated  by  setting  the  time 
constant  of  the  guidance  system,  r2»  to  .1  seconds  in  the 
truth  model.  t2  in  the  filter  was  set  to  0.3  seconds.  The 
fourth  order  filter  was  used  and  only  the  dynamic  states  of 
the  missile  model  were  estimated. 
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Fig.  117.  RANGE  SI  OMRS  FOURTH  ORDER 
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Sens itivltv  Analysis  ( M = 2.) 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 
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This  set  of  plots  was  generated  by  setting  the  mass  in 

the  truth  model  to  2 slugs.  The  fourth  order  filter  was 

used  with  the  mass  of  the  missile  in  its  model  set  at  4. 

Both,  filter  and  truth  model  used  an  S (for  M/S)  of 
2 

.137  ft  . Only  the  dynamic  states  of  the  missile  model  were 
estimated. 
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Fig.  123.  X-VELOCITY  SIDHRS  FOURTH  ORDER 
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Sensitivity  Analysis  (M  = 8. ) 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 
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This  set  of  plots  v/as  generated  by  setting  the  mass  in 

the  truth  model  to  8 slugs.  The  fourth  order  filter  was 

used  with  the  mass  of  the  missile  in  its  model  set  at  4. 

Both,  filter  and  truth  model  used  an  S (for  M/S)  of 
? 

.137  ft  . Only  the  dynamic  states  of  the  missile  model  were 


estimated. 
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Fig.  135.  ANCLE  MEASUREMENT  S1CMAS  FOURTH  ORDER 
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Fig.  137.  RANGE  SIGMAS  FOURTH  ORDER 
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Fig.  138.  RANGE  RATE  FOURTH  ORDER  MISSILE 
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Fig.  I4l.  LRT  RCCELERRTION  SIDMRS  FOURTH  ORDER 
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First  Order  Missile  Filter  (r^  set  equal  to  .85) 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 

V1  (0)  = 1225.?  fps 
mx 

^(0)  = ^ . 363345  radians 

R( 0)  = 10000.  feet 

R( 0)  = -2122.  fps 

a^O)  = 0. 

n(0)  = 4.5 

T^( 0)  = .85  seconds 

M/S(0)  =29.197  slugs/ft2 


These  plots  were  generated  by  the  first  order  filter 
with  a time  constant  of  .85  seconds.  The  value  of  was 
found  from  an  iterative  process  by  comparing  filter  per- 
formance for  various  values  of  t^..  .85  seconds  was  found 

to  produce  the  least  error  in  the  dynamic  states  of  the  mis- 
sile model.  Only  the  dynamic  states  were  estimated  in  this 
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Fig.  143.  X-VELOCITY  SICMflS  FIRST  ORDER 
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Fig.  144.  RNOLE  MERSUREMENT  FIRST  ORDER  MISSILE 
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Fig.  145.  ANGLE  MEASUREMENT  SIGMAS  FIRST  ORDER 
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Fig.  14?.  RANGE  SIGNRS  FIRST  ORDER 


HEflN  ERRO* 


t£ 

hJ 

UJ 

ID 

-J 

*-• 

*— « 

-J 

to 

*-* 

O 

in 

UJ 

u_ 

*-* 

%- 

S. 

in 

c 

ac 

03 

K 

X 

W 

UJ 

UJ 

Ui 

V- 

o 

ti 

3 

a 

<*: 

z 

o£ 

o 

cc 

V- 

i— 

a; 

Li 

V— 

a 

W 

s 

r> 

*■* 

a 

u. 

w 

»~ 

1 ( 

CO '21  00‘6 

03S 


Fig.  149.  RRNOE 


in 

O) 

in 


RTE  S I OMRS  FIRST  0R0ER 

12 


riESN  ERRIM 
LflT  ACCELERATION 


00 ’9 

00  ’£  OO'O  00*c- 

00*9° 

aaniiNouw  uouio 

Fig.  150.  LOT  ACCELERATION  FIRST  OROER  MISSILE 


183 


5.95 


;lflNCE 

/FILTE 


150. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

l.E- 

6 0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

500. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

200. 

0. 

0. 

0. 

0. 

k = 

0. 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

.01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

_ 0* 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

j 

This  set  of 

plots  was 

generated 

by  the 

first 

order  fil 

ter  when  estimating  n,  which  was 

initialized 

at 

3 

in  the 

filter. 

. The 

true 

value  of 

n in 

the 

truth  model 

was  4.5. 

The  other  parameters  in  the  filter  model  were  not  estimated. 
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This  set  of  plots  was 

generated  by  the 

first 

order  fil 

ter 

when  estimating 

n,  which  was 

initialized 

at 

6 

in  the 

filter. 

The 

true  value  of 

n in 

the 

truth  model 

was  4.5. 

The  other  parameters  in  the  filter  model  were  not  estimated. 
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Fig.  165.  X-VELOCITY  FIRST  ORDER  MISSILE 
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Fig.  166.  X-VELOCITY  SICNRS  FIRST  ORDER 
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Fig.  170.  RRNBE  FIRST  OROER  MISSILE 
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Fig.  171.  RANGE  SIGMAS  FIRST  ORDER 
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Fig.  172.  RBNOE  RRTE  FIRST  ORDER  MISSILE 
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Fig.  174.  LBT  BCCELERRTION  FIRST  ORDER  MISSILE 
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Fig.  177.  PRO-NRV  CONSTANT  FIRST  OROER  MISSILE 
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Xf  Estimation  - initialized  at  1 . 5 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 

1225-7  fps 
4.363345  radians 
10000.  feet 
-2122.  fps 
0. 

4.5 

1.5  seconds 
29.197  slugs/ft2 


= 

0/0)  = 
R(0)  = 

R(o)  = 
a^O)  = 
n(0)  = 

Tf(0)  = 

m/s( 0)  = 


These  plots  were  generated  by  estimating  with  the 
first  order  filter,  t^.  was  initialized  at  1.5  seconds 
with  its  truth  model  value  defined  as  0.85  seconds  The 
truth  model  value  was  determined  from  an  interative  search 
as  described  in  Chapter  IV. 
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Fig.  179.  X-VEIOCITY  FIRST  QROER  MISSILE 
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Fig.  181.  RNDLE  MERSUREMENT  FIRST  ORDER  MISSILE 
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Fig.  I83.  RRNOE  FIRST  ORDER  MISSILE 
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Fig.  185.  RANGE  RATE  FIRST  ORDER  MISSILE 
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Fig.  188.  LRT  RDCELERBTIQN  S10MRS  FIRST  ORDER 


22? 


Xf  Estimation  - Initialized  at  . 3 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 


V;V(0) 

mx 

= 1225.7  fps 

e^o) 

= 4.363345  ra< 

R(o) 

= 10000.  feet 

R(o) 

= -2122.  fps 

^(o) 

= 0. 

n(0) 

= 4.5 

Tf(0) 

= .3  seconds 

M/S ( 0 ) = 29.197  slugs/ft" 


r.E-5  o.  o. 

R = 0.  500.  0. 

0.  0.  100. 


100.  0.  0. 

0.  l.E-8  0. 

0.  0.  101 
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SL  = 
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0.  200.  0. 

0.  0.  .1 
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0. 
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0. 
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0.  0.  0.  0. 


0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  .001  0. 

0.  0.  0. 


These  plots  were  generated  by  estimating  with  the 
first  order  filter,  was  initialized  at  0.3  seconds  in 

the  filter  with  its  truth  model  value  defined  as  O.85 
seconds.  The  truth  model  value  was  determined  from  an 
iterative  search  as  described  in  Chapter  IV. 
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Fig.  191.  X-VELDCITY  FIRST  ORDER  MISSILE 
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Fig.  192.  X-VE10C1TY  S1DMRS  FIRST  ORDER 
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Fig.  193.  ANDIE  MEASUREMENT  FIRST  ORDER  MISSILE 
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Fig.  195-  RANGE  FIRST  ORDER  MISSILE 
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Fig.  196.  RANGE  SIGMAS  FIRST  ORDER 
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Fig.  197.  RANGE  RATE  FIRST  ORDER  MISSILE 
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Fig.  201.  TflU  PARAMETER  FIRST  ORDER  MISSILE 
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M/S  Estimation  - M/S  Initialized  at  45. 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 

VL(0)  = 1225  * 7 fps 

9(0)  = 4.363345  radians 

T 

R( 0)  = 10000.  feet 

R(0)  = -2122.  fps 

aL(0)  = 0. 

n(0)  = ^.5 

T^.(0)  = .85  seconds 

M/S( 0)  = 45.  slugs/ft2 
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These  plots  were  generated  by  the  first  order  filter 
when  estimating  the  M/S  ratio.  The  initial  value  of  this 

p 

ratio  set  in  the  filter  was  45.0  slugs/ft  . The  true  value 

p 

was  set  at  29.197  slugs/ft  . This  parameter  along  with 
the  five  dynamic  states  of  the  missile  model  were  the  only 
estimates  of  the  first  order  filter  in  this  case. 
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Fig.  203.  X-VELOCITY  FIRST  ORDER  MISSILE 
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Fig.  204.  X-VELOCITY  SIGMAS  FIRST  ORDER 
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Fig.  206.  RUBLE  MEASUREMENT  FIRST  ORDER  MISSILE 
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Fig.  203.  RfiNDE  FIRST  ORDER  MISSILE 
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Fig.  210.  RRNBE  RATE  FIRST  ORDER  MISSILE 
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Fig.  212 . LRT  ACCELERfiTlON  FIRST  ORDER  MISSILE 
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Fig.  213.  LfiT  RCCELERRTIQN  SIDMflS  FIRST  ORDER 
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Fig.  214.  VRRIRNDE  CONVERGENCE 
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Fig.  216.  MfiSS/SURF  SIGMRS  FIRST  ORDER 
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M/s  Estimation  - M/S  Initialized  at  15- 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 


yL‘0)  = 
0(0)  = 

R(0)  = 

R(0)  = 

«l(0)  = 

n(0)  = 

Tf(0)  = 

M/S(0)  = 


1225-7  fps 
4.363345  radians 
10000.  feet 
-2122.  fps 
0. 

^•5 

. 85  seconds 
15  slugs/ft2 
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0. 
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5. 

260 


These  plots  were  generated  by  the  first  order  filter 
when  estimating  the  M/S  ratio.  The  initial  value  of  this 

p 

ratio  set  in  the  filter  was  15.O  slugs/ft  . The  true  value 

p 

was  set  at  29-197  slugs/ft  . This  parameter  along  with  the 
five  dynamic  states  of  the  missile  model  were  the  only  esti- 
mates of  the  first  order  filter  in  this  case. 
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Fig.  217.  X-VELOCITY  FIRST  ORDER  MISSILE 
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Fig.  220.  RNGLE  MEASUREMENT  FIRST  ORDER  MISSILE 
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Fig.  221.  ANCLE  MEASUREMENT  SIGMAS  FIRST  ORDER 
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Fig.  222.  RANK  FIRST  ORDER  H1SS11E 
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Fig.  225.  RANGE  RATE  SIGMAS  FIRST  ORDER 
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Fig.  226.  LAT  ACCELERATION  FIRST  CROER  N1SS1LE 
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Fig.  227.  LRT  ACCELERATION  S10MRS  FIRST  ORDER 
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Fig.  229.  MRSS/SURF  FIRST  OROER  MISSILE 
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The  initial  state  estimates  and  the  tuning  parameters 


for  this  case  are 


V*  (0) 
mx 

e/o) 


R(o) 

R(o) 


«l(°) 

n(0) 

Tf(0) 

m/s( 0) 


1225.7  fps 
4.3633^5  radians 
10000.  feet 
-2122.  fps 
0. 

6. 

.3  seconds 
29.197  slugs/ft2 
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001 
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These  results 

were 

produced 

when 

estimating 

n and  t^. 

l 

together  with  the 

five 

dynamic  states. 

Since 

these  para- 

me ters 

only 

appeared  together  in 

the  model  of 

the 

guidanc 

strip,  they  were  considered  the  most  difficult  combination 
to  estimate  simultaneously,  n was  initialized  at  6. , and 
Tf  was  initialized  at  0.3  seconds.  The  true  values  of  n 
and  Tf  were  4.5  and  O.85  seconds  respectively. 


MEAN  EIWOR 


Fig.  231.  X-VELOCITY  FIRST  OROER  MISSILE 
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Fig.  235-  RRNOE  FIRST  ORDER  MISSILE 
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Fig.  242.  PRO-NRV  CONSTANT  S1GMRS  FIRST  ORDER 
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Fig.  243.  TfiU  PARAMETER  FIRST  ORDER  MISSILE 
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Fig.  244.  TAU  PARAMETER  SIDMAS  FIRST  ORDER 
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n and  Estimation  (low-g  scenario) 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 

V*  ( 0)  = 1225.7  fps 
mx 

9t(0)  = ^.363345  radians 

R( 0)  = 10000.  feet 

R(0)  = -2122.  fps 

a1(0)  = 0. 

n(0)  = 3- 

Tf(0)  = 1.5  seconds 

M/S( 0)  = 29.197  Slugs/ft2 
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These  plots  were  generated  as  a test  case  to  determine 
if  the  parameters  could  be  estimated  for  a low-g  scenario. 
Only  n and  were  estimated  along  with  the  five  dynamic 
states  of  the  missile  model.  These  plots  can  be  compared 
with  Figures  231  through  244  which  were  generated  using 
the  high-g  scenario.  The  n and  Tf  in  the  filter,  for  this 
set,  were  initialized  at  3.0  and  1.5  seconds  respectively. 
The  true  n was  4.5  and  thr  true  value  for  was  defined  as 
0.85  seconds. 


293 


sro  so  *0 

00 T0  80*0- 

sro-° 

,.01  * 

SNuiauy 

Fig.  247.  ANGLE  MEASUREMENT  FIRST  ORDER  MISSILE 
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Fig.  249.  RANGE  FIRST  ORDER  N1SSILE 
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Fig.  254.  LRT  ACCELERATION  SIGMAS  FIRST  QROER 
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Fig.  256.  -RO-NRV  CONSTANT  SICMAS  FIRST  ORDER 
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n,  jf,  and  M/S  Estimation 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 


yL(o> 

0(0) 

R(0) 

R(0) 


1225.7  fps 
4.3633^5  radians 
10000.  feet 
-2122.  fps 


£^(0)  = 0. 

n(0)  = 6. 

Tf(0)  = .3  seconds 

M/S(0)  = 15.  slugs/ft2 
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These  plots  were  generated  to  demonstrate  the  first 
order  filter's  ability  to  estimate  all  three  parameters  and 
all  five  dynamic  states  simultaneously.  The  states  were 
not  initialized  with  any  error.  The  parameter  were  ini- 
tialized in  the  filter  as  follows 

n(0)  = 6.0  , true  value  = ^.5 

TfCO)  =0.3  seconds  , true  value  = .8 5 seconds 

M/S(0)  = 15.  slugs/ft2  , true  value  = 29.197  slugs/ft2 
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Fig.  260.  X-VEL0C2TY  SIDKflS  FIRST  ORDER 
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Fig.  262.  RNOLE  MERSUREMENT  FIRST  ORDER  HISSILE 
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Fig.  263.  ANGLE  MEASUREMENT  SIGMAS  FIRST  ORDER 


314 


00 ’0£ 

oo*st 

00 -0 

oo*si- 

00 ‘0E-0 

1333 

Fig.  264.  RANGE  FIRST  ORDER  MISSILE 
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Fig.  268.  LflT  ACCELERATION  FIRST  ORDER  MISSILE 
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Fig.  273.  7RU  PRRRMETER  FIRST  ORDER  MISSILE 
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Fig.  2 74.  TAU  PARAMETER  SIGMAS  FIRST  ORDER 
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Fig.  275.  MRSS/SURF  FIRST  OROER  MISSILE 
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Fig.  276.  MRSS/SURF  SIGMAS  FIRST  ORDER 
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n,  Xf  M/S  Estimation  with  Initial  State  Errors 

The  initial  state  estimates  and  the  tuning  parameters 
for  this  case  are 


aL(0)  = 15.  g's 

n(0)  =6. 


Tf(0)  = .3  seconds 

M/S( 0)  =15.  slugs/ft2 


This  was  the  final  evaluation  of  the  first  order  filter. 
The  states  and  parameters  were  both  initialized  with  large 
errors.  The  initial  values  of  the  filter's  states  are 
given  on  the  previous  page.  The  true  values  for  each  are 
listed  below. 


VL(0)  = 1225-6  fps 

^(0)  = 4.3633^5  radians 

R(0)  = 10000.  feet 

R(0)  = -2122.  fps 

aL(0)  = 0. 

n(0)  = 4.5 

Tf(0)  = .85  seconds 

M/S(0)  = 29.197  slugs/ft2 
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Fig.  2?8.  X-VEL0C1TY  SIOMRS  FIRST  ORDER 
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Fig.  280.  ANCLE  MEASUREMENT  FIRST  ORDER  MISSILE 
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Fig.  282.  RRNOE  FIRST  OROER  MISSILE 
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Fig.  283.  RRNBE  81 OMRS  FIRST  ORDER 
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Fig.  284.  RANGE  RATE  FIRST  OROER  MISSILE 
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Fig.  285.  RANGE  RRTE  S1GHR3  FIRST  ORDER 
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VRRIRNCE 


Fig.  288.  VARIANCE  CONVERGENCE 
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Fig.  289.  PRO-NAV  CONSTANT  FIRST  ORDER  MISSILE 
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Fig.  290.  PRO-NAV  CONSTANT  SIBMRS  FIRST  ORDER 

343 


08*0 


00*lP 


09'0  0t*0  02*0 

30niINOUU  V0VU3 

Fig.  291.  TRU  PRRRMETER  FIRST  ORDER  MISSILE 
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Fig.  292.  TRU  PARAMETER  SIGMAS  FIRST  ORDER 

345 


00 't 


00'C  00  'Z  OO’I 

30fli INOUW 


OO’lP 


